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Abstract. We study the Dirichlet problem associated to the equa- 
tion for self-similar surfaces for graphs over the Euclidean plane 
with a disk removed. We show the existence of a solution provided 
the boundary conditions on the boundary circle are small enough 
and satisfy some symmetries. This is the second step towards the 
construction of new examples of complete embedded self similar 
surfaces under mean curvature flow. 



1. Introduction 

This paper is the second one of a series of three articles describing 
the construction of new examples of complete embedded self-similar 
surfaces under mean curvature flow [8][7j. Our general strategy is in- 
spired by Kapouleas' article [3] and is described in detail in the previous 
installment. Let us recall it briefly: to construct a new self-similar sur- 
face, we take two known examples suitably positioned and replace a 
neighborhood of their intersection with an appropriately bent scaled 
Scherk's singly periodic surface. The procedure is called desingulariza- 
tion. 

The resulting surface is not smooth; however it is a good approximate 
solution. The next task is to find a small function whose graph over it 
satisfies the self-similar surface equation 

(1) H + X-u = 0, 

where H is the mean curvature and v is the normal vector so that the 
mean curvature vector is H = Hv. The sign of H is chosen so that 
the mean curvature of a convex surface is positive. Before consider- 
ing graphs of functions on the entire surface, we have to work locally 
and study the Dirichlet problems with small boundary conditions cor- 
responding to ([T]) for graphs of functions on the different pieces. 
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In the first part [8], we found that a perturbation of an appropri- 
ately scaled and bent Scherk surface satisfies (jTJ and can be used to 
desingularize the intersection of a cylinder and a plane or a sphere and 
a plane. In this article, we study the Dirichlet problem corresponding 
to (JT]) for functions on the outer plane, which is the Euclidean plane 
with a disk removed. We show that the Dirichlet problem possesses 
a solution if the boundary conditions on the circle are small enough. 
Our methods here are different from the one used in Part I. Since the 
equation ([T]) can be written with a single set of coordinates on the outer 
plane, our tools come mainly from partial differential equation theory, 
with some interior estimates established by Ecker and Huisken [TJ for 
the Mean Curvature Flow. The third article is to discuss the gluing of 
the solutions on the different pieces in a manner to obtain a smooth 
complete embedded self-similar surface. 

1.1. Main result. Let 2 > R > \/3/2 and N > 5, consider a function 
/ E C 4 ([0,27r]) satisfying the symmetries 

(2) f{ e) = -f(-6) = f(n/N-6). 

Theorem 1. There is an £q > depending on R and N such that for 
any f : [0,2n] — > R with ||/||c 4 [o,27r] = £ < £ o an d f satisfying the 
symmetries above, there is a function u on Q = R 2 \ Br such that 

(3) the graph of u over Q satisfies H + X ■ v = 

(4) u = f on dB R , 

(5) u(r, 6) = -u{r, -9) = u{r, n/N — 6) for r > R,8 E [0, 2vr) 

In addition, We can choose the constant Eq uniformly for all R such 
that 2 > R > y/3/2. 

Here, u = f on OBr means that u(R, 6) = f(8) in polar coordinates. 
This is a slight abuse of notation which does not induce confusion. Note 
that the conditions ([5]) imply that u is odd with respect to rotations 
of 180 degrees with respect to the x-axis and even with respect to 
reflections across the planes 9 = tt/2N + kn/N, which are exactly the 
symmetries imposed in Part I [5]. 

Computing the mean curvature and the normal vector v of the graph 
of u using coordinates in R 2 , we find that ([3]) is equivalent to 

(6) £(u)=g ij (Du(Z))D ijU (t;)-Z-Du(Z) + u(Z) = 0, £ E tl, 

where g^(Du) = 5^ — 1H !" D „p ■ This is our main equation. Note that our 
domain is unbounded and that the coefficient in front of u is positive, 
so the standard maximum principle is not applicable directly. 
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We first study the properties of the linear operator C associated to 
£. From solutions to the linear problem, we then construct a sub and 
supersolution to the quasilinear equation £u = 0. 

The solution to the elliptic equation ([6]) is found as a limit for time 
r going to infinity of a solution u to the parabolic equation d T u = 
£u. For initial conditions with bounded gradient, standard theory on 
parabolic equations assures the existence of a solution to the initial 
value problem on a short time interval [0, To). The constructed sub and 
supersolution serve as barriers for solutions of the parabolic equation 
d T u = £u and allow us to control the growth of u at infinity in the 
space variable for all time. Exploiting the close relation between the 
parabolic equation d T u = £u and the mean curvature flow, we use 
interior estimates from the latter to bound derivatives of our solution 
to the former away from the boundary. Loosely speaking, this gives 
us control of the derivatives outside of a bounded annulus around dfl. 
To bridge the gap, we invoke standard parabolic theory on bounded 
sets and obtain uniform estimates for derivatives of u for all r £ [0, r ). 
This implies that the solution has to exists for all time r £ [0, oo). We 
conclude that a subsequence of «(■, r) tends to a solution of £u = as 
r goes to infinity by proving a monotonicity formula. 

Some sections of this paper appeared in the author's thesis, written 
under the direction of Sigurd Angenent at the University of Wisconsin- 
Madison. The author is indebted to him for invaluable discussion. 



2.1. Mean curvature flow. The properties of solutions to ([3]) are 
closely related to properties of solutions to the mean curvature flow 
(MCF), which will be used intensively. Let us therefore define it. 

Let X(-,t) : M 2 — > R 3 be a one parameter family of immersions of 2- 
dimensional smooth hypersurfaces in R 3 . We say that M t = X(M 2 ,t) 
is a solution to the mean curvature flow if 



is satisfied for some initial data X . Here H(p, t) is the mean curvature 
vector of the hypersurface M t at X(p,t). 

If the hypersurfaces M t can be written locally as graphs of a function 
v(-,t) over a domain in the xy- plane, the quasilinear equation 



2. Definitions 



(7) 
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X(p,t) = H(p,i), peM,t>0 
X(p,0) = X (p) 
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is equivalent to (J7j) up to tangential diffeomorphisms. 

2.2. Weighted Hilbert Sobolev spaces H, V and Vo. We consider 
the Hilbert space 

H = L 2 (Q,dm(0), 
where Q is an open subset of R n (not necessarily bounded) and m is 
the Gaussian measure dm{^) = e~ 2 1^1 d£. We write W/Wh and (f,g)n 
for the norm and inner product of H. We also define the Hilbert space 

V = {f e H : dif 6 H for i = 1, . . . , n}. 

The inner product of V is (f,g)v = J n {j 9 + D f ■ Dg}dm(£ i ) . Let Vo be 
the closure of C^°(Q) in V. Even when Q is unbounded the inclusion 
Vo C H is compact. 

3. The Linear Operator 
Let us define sections fl Rt N of the outer plane by 
(9) n RiN = {(r,6) : R<r, -tt/N < 9 < ir/N} 

and their corresponding Hilbert spaces 

V 0;Sym (^R,N) = {ue V (n R , N ) \ u(r, 9) = -u(r, -9) = u(r, tt/N - 9)} 
H Sym (n R , N ) = {ue H(n R>N ) I u(r, 9) = -u(r, -9) = u(r, vr/iV - 9)} 

In this section, we study the properties of C, the linear operator asso- 
ciated to £ defined in (jSJ), 

C(u) = Ait — £ • Du 4- u. 

3.1. Eigenvalues and eigenfunctions. Consider the operator 

= Am — £ • Du — u 

on f2 ,AT. The compactness of the inclusion Vo(fio,iv) C H{VLq,n) yields 
the following theorem: 

Theorem 2. The operators £_ and £ have a countable set of eigen- 
values having no limit point except possibly A = oo and a corresponding 
H -orthogonal basis of eigenfunctions tpt G Vo(^o,jv)- 

Proof. Theory on compact bounded self-adjoint operators in Hilbert 
spaces gives us the result for The result for C then follows imme- 
diately. □ 
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Consider eigenf unctions of the form f{r)g{9) satisfying the symme- 
tries 

(10) g(e) = -g(-0)=gWN-O). 

There is an eigenvalue A G R such that 

1 1 

£{f{r)g{0)) = f rr g + -f r 9 + -^fgee - rf r g + fg = Xfg. 

Separation of variables and the symmetry conditions (flOl) on g imply 
that the only possibilities for g are g = C sin(m#) with m = kN, k G N. 
Therefore 

gee = ~m 2 g 

C (f) = frr + ~fr ~ rf r + (1 - A)/ - f = 0. 

We have 

C (r n ) = {n 2 - m 2 )r n - 2 + (1 - A - n)r n . 

jCq is a linear map on the vector space V = {r m F(r 2 ) \ F polynomial } 
and it is upper diagonal in the basis {r m , r m+2 , r m+4 , . . .}. The operator 
jCq has a zero eigenvalue if one of the entry on the diagonal vanishes, 
i.e. if 1 — A — (m + 21) =0. We therefore found some of the eigenvalues 
of our original operator C, 

(11) A = 1 - (kN + 21), kN,(6 {0} U N, 

where N = {1,2,...} is the set of natural numbers. For each value of 
k, I the corresponding eigenfunction is of the form 

(12) r kN P Kl {r 2 ) sm{kN6) 
where Pk,i is a polynomial of degree I. 

We now show that the eigenvalues values given by ( TTTj) are in fact 
all the possible eigenvalues of C. This is done by proving the following 
theorem: 

Theorem 3. The set S of eigenf unctions (1121) corresponding to the 
eigenvalues (HID forms a basis for the Hilbert space Hs ym {Qo,N) ■ 

Before we start the proof, note that for each fixed fcsN and s/l, 

(r kN P k jsm(kN6),r kN P Ks sm(kNe)) H{nQN) = 
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since they are eigenf unctions corresponding to different eigenvalues. 
Hence, J^ Q Pk,i(r 2 )Pk,s( r2 ) r2kN &~ r ^ r dr = 0. By the change of vari- 
able p = r 2 , 

POO 

(13) / P k , l (p)P k , s (p)p kN e-^ 2 dp = 0. 

Jp=0 

The polynomials P^i, which are closely related to the general Laguerre 
polynomials, are therefore orthogonal in L 2 (0, oo; p kN e~ p ^ 2 dp), the L 2 - 
space weighted by p kA r e~~ p ' 2 dp. In addition, they span the functions 
p n , n G N. The fact that {Pk,i}i^{o}uN is a maximal orthogonal set in 
L 2 (0, oo; p kN e~ p / 2 dp) follows from the theorem p. 333 in [9] which is 
recalled below. 

Theorem 4 (Theorem p. 333 [9J). The orthonormal sequence of poly- 
nomials attached to a mass distribution function p(x) on the interval 
(a,b), is complete in L 2 (a,b;p) whenever there exists a number r > 
such that the integral 




e rlxl dp{x) 



exists. 

Proof of Theorem 0. It suffices to show that the set of all finite linear 
combinations of members of S is dense. Let f(r,9) G Hs ym (Qo,N)- 
Since / is periodic in 9 and f(r, ■) G L 2 (—tt/N,tt/N) for almost every 
r, we can expand f(r, •) in Fourier series for almost every r G (0, oo). 
The symmetries of / imply f(r,9) = X^gN a j( r ) siix(jN0). Using the 
change of variable p = r 2 , we write 

(14) f(p, 9) = f(^p, 9) = J2 HP) MJN9) for a.e p G (0, oo). 
The if-norm of / is finite, therefore 

p poo ptt/N 

/ \f(r,9)\ 2 e-W 2/2 d£= / \f(r,9)\ 2 d9 e~ r2/2 r dr 

JQq^n Jr=0 J8=-tt/N 

= W £W C ~^=^E/ b 2 (p)e^ 2 dp <oo. 

Hence, the functions 6,'s are in L 2 (0, oo; e~ p/2 dp) for j G N. For every 
e > 0, there is a J(e) G N so that 

POO 
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We can approximate bjp~i N / 2 G L 2 (0, oo; fP N e~ p l 2 dp) by a linear com- 
bination of Pj t i, since {-Pj,z}«e{o}uN is a maximal orthogonal set . More 
precisely, for every e > and every j, there is a linear combination of 
Pj/s denoted by Qj )E so that 

Define Q £ to be the linear combination of elements of S given by Q £ = 
Eo<i<j(e)Q^ N/2 sHjN6)- Then 

11/ - Qe\\H(n 0iN ) 

• 7p -° \i>J(e) 0<j<J( £ ) / 

< e/2 + 2 ~ J ~ l£ < £ - 

0<j<J(s) 

This shows that every function in Hs ym (fio t N) can be approximated 
arbitrarily closely by a linear combination of functions of S. □ 

3.2. Poincare inequality. As an immediate consequence of Theorem 
[3l we have the following inequality for u G C£°(flo,.iv): 

= -(Au-£- Du-u,u) H > -A ||u||h = (-W + l)Nlzr, 

where A is the lowest eigenvalue of the operator — This implies a 
Poincare inequality 

(15) \\u\\v(cv) > VN + 1 ||«||H(no, « e Vb(fi') 
for every domain Q' C f^o.jv- 

3.3. Dirichlet problem. Let fi^jv be defined as in (jHJ). 
Lemma 5. For any g G H(Q^n), the equation 

(16) £(u) = Au-£- Du + u = g 

possesses a weak solution u G Vo(f2R,jv)- 

Proof. It is a classical variational argument for the functional J(u) = 

L R , N {h\ D <0\ 2 - HO 2 + <M0> over « G K,(^ )JV ) and 

using the Poincare inequality (jT5l) . □ 
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3.4. Maximum Principle. Let Q' be any domain that is a (not nec- 
essarily proper) subset of fi^jv- 

Definition 6. We say that u G V(fi') satisfies u < on dfl' if its 
positive part u + = max(w, 0) G Vq(Q'). 

Theorem 7 (Maximum Principle). Let N > 5. Suppose u G V(O') 
satisfies 

£w = Am — £ • Dm + u > m fi' 
u < on dVt' , 

then u < in Q' . 

Proof. First note that ^ := rcos 8 is a solution of £^ = 0. By the 
hypotheses, we have 

/ Du ■ Dv - uv dm < 0, for v > 0, v G C^tt'). 
Jw 

Define the function ip such that u = if) (p. Since l/if> and l-D^I are 
uniformly bounded, </? G V. Moreover, ip satisfies 

/ ipDip ■ Dv — vDip ■ Dip + Dip ■ D(ipv) — ifupv dm < 0. 
Jn> 

Since Cip = and <pv G Vo, we have J n , Dip ■ D(<pv) — ifxpv dm = 0. 
Hence, 

(17) / ipD<p- Dv -vD<p- Dip dm<0. 

Jn> 

Take v = max(<p,0) in (|17|) . Note that we can do this since u + G Vo 
and ^ 7^ 0, so v G Vo- Therefore 

ipDv ■ Dv dm < I vDv ■ Dip dm 
Jn> 

since Dv = Dip on the support of v. Using i?cos ^ < ip and \Dip\ < 1, 
we get 



7T , ,„ , 



and 



^ COS ~2N dm < I \Dv\v dm, 



7T 2 I, 

-Rcos— < HhIpIIh 



by the Cauchy Schwarz inequality. Our Poincare inequality ( TT5|) implies 
that 

7T 1 

Rc0S 2N W Dv Wh < P u ll^-7=ll^lli?- 
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A simple computation shows that for iV > 5, R cos ^ — -^ > ^75 cos 7^— 
h= > 0. Therefore 

\\Dv\\ H = Q. 

But v G Vo(Q'), so v = 0; this means that <^ < and w < in fi'. □ 

3.5. Boundary value problem. Let F G C 4 ([0,27r]) with the sym- 
metries (ED : F(9) = -F{-6) = F(n/N -8). We will use the notation 
Kq = ||i 7 '||c 4 ([o,27r]) throughout Section [21 

Definition 8. We say that u G V(Q) (V(Qr,n)) is a solution to the 
problem 

Cu = 0, in Q (fl,R t N resp), u = F on dfl (OQr^ resp. ) 
if u satisfies 

Cu = 0, in Q {Qr,n resp.), u — if) G Vo(Q) (Vo(Qr,n) resp.), 

where if) is a function in C 4 (f2) (C^^r^) resp.) with ip(R,9) = F{9) 
and if>{r,6) = -if)(r,-9) = if;(r, tt/N '- 9) for r > R,9 G [0,2n), 
(r>R,9e [-tt/N, n/N] resp.). 

It follows from Lemma [5] and the Maximum Principle (Theorem [7j) 
that the solution u to Cu = in VIr^n, u = F on dflR^ exists and is 
well defined, in other words it is unique and independent of the choice 

of if). 

Denote by U\ the extension of u to Q given by 
U 1 (r,0)=u(r,9) mn R>N 
U x (r, 9) = -U x (r, -9) = U x (r, tt/N - 8) in Q. 
The function U\ is a solution to boundary value problem 

(18) CUi = AC/i - £ ■ DU X + U x = in fi, C7i = F on dVt. 
Moreover, it is smooth by standard elliptic theory. 

3.6. Behavior of solution at infinity. Let u be the solution to the 
problem Cu = in Qr^n, u = F on dQ,R t N. We can choose if) in 
Definition [S] satisfying Halloo < Halloo with if) = in Qr^n \ B R+2 . Let 
a = 3 1| -P || oo and consider the function 

W\ = u — ar cos 9. 

Since N > 5, a = > Jf^fe) therefore wi < on OQr^. 

It follows from the Maximum Principle that u>i < in £Ir } n, hence 
u(r, 9) < ar cos 9 in O^at- An similar argument with w 2 = —ar cos9 — u 
gives a lower bound on u{r,9), so 

(19) \u{r,9)\ <arcos{9) <3K r on f^Rjv- 
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Therefore, we also have 

(20) 1^(01 <3tf |£|, (efi. 

3.7. Bounds on the first and second derivatives of the solution. 

Let Ui be the solution to ffl8|) from section 13.51 

Lemma 9. There is a constant K\ independent of U\ and F so that 

(21) {DUxl + ltWD 2 ^ <if 1 || J P|| £7 *[o,jhr] 

Proof. Step 1: Estimate away from the boundary. 

The function U\ is smooth in Q by elliptic theory. A computation 

shows that the function 

x 



)(x,t) = y/l-t Ui 



satisfies the heat equation v t = Av. The theory for parabolic equations 
gives us the following estimates on the derivatives of v, 

(22) |I/^( P) 0)|<%) sup \v(x,t)\, 

« B s (p)x(-s 2 ,0) 

where C{m) is a constant independent of v or s and as long as v exists 
in B s (p) x (— s 2 ,0). This is true provided 

minseB.fr) \x\ = \p\ - s 
max (6( _ s2 , 0) y/2(l -t) ^2(1 + s 2 ) 

In particular, it is true if \p\ > ^pi? and s = |p|/10. From (1221) . we get 
\Dv{p,0)\ < y sup \v(x,t)\. 

\P\ B s (p)x(- S 2,0) 

The last factor of the right hand side can be estimated using equation 
(j20j) to obtain 

sup \v{x,t)\ < V2 K (\p\ + s). 

B e (p)x(-s*,0) 

Using the notation £ = pj \/2, we have 

|-Dt/x(0l = V^pu(p,0)| < 22 C(1)AT for \i\ > 5R/3. 
Similarly, we get an estimate on the second derivative 

|D 2 fMOI = 2\D 2 v( Pl 0)\ < 22 °3 2) Ko for |£| > 5i?/3. 

Therefore 

(23) \DU l \ + |£||£> 2 £/i| < C||F|| L oo for |f | > 5/2/3. 
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Since R < 2, the estimate ( |23|) is valid for |£| > 4 in particular. 

Step 2: Estimates up to the boundary. 
Consider the equation CU\ = in the annulus A Rji = £> 4 \ B R . The 
domain A R;i is bounded so global regularity results for elliptic equations 
from [3] imply that 

(24) ||t/i|k« ( A fl , 4 ) < CfllE/Ji*^) + ||V||^(^ 4 )) S 

where ip is a function in W^'^Ar^) such that ^ — U\ G Wo' 2 (Ar,4)- 

Note that we can bound the first four derivatives of U\ for |£| = 4 
by an argument similar to the one in step 1, 

< C\\F\\ L oo for i = 1,...,4 and £ G d5 4 - 

It is then clear that we can find a constant C independent of F and 
£7i, and a function -0 for which 

1p\dB R = F, lp\dB4, = UxlsBt, 

(25) U\\c*(A R , 4 )<C\\F\\ C 4 [0M . 
With this choice of ip, equation (1241 gives us 

(26) ll^ilk^^) < L7(||^|| W 4,2 (Afl4) + WU-lW^a^). 
Using Sobolev's inequality, ( 1251) . ( 1261) and ( 1201) . we get 

(27) ||C/i|| C2 (A H , 4 )<C||F|| C 4 [0)27r] . 

Combining (1231) and (1271) . we obtain the desired result. □ 

4. Finding sub and supersolutions 

This section is devoted to the construction of a subsolution w_ and 
a supersolution u + to the problem ©-([4]), or equivalently, to 

DjuD~uDjjU 

(28) £u = Au-£- Du + u- 3 3 



1 + \Du\ 2 
u\ d n = sF 

with small e. Let us discuss the strategy for finding a supersolution 
first. Since the boundary data is of order e, we write u = elf, 

CU ~ £ l+e*\DU\* ' Ulm ~ F 
and decompose U into the two terms U — U\ + £ 2 ?72, where 

(29) £fA=0, C/ilan^F, 

(30) IT, . tf a | m = 0. 
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We discussed the existence of such a function U\ in section [331 Choos- 
ing an appropriate function C/ 2 is done below. 

4.1. Preliminary computations for finding a supersolution. In 

order to get a supersolution, we want U2 to satisfy CU2 — ^j^r^m~ — 
0. Roughly estimating U by U\ and bounding the derivatives of U\ 
using Lemma [9], we are looking for a C/j so that 

where Ci = i^iHi^Hc 4 ^^- Define v a = a (r — \ . A simple compu- 
tation gives 

aR 2 a 



Cv a = + -(1 - 2iT), ?;„ = on 



This makes t> a a perfect candidate for f/ 2 - 

Remark 10. The term a(l — 2R 2 )/r has to counterbalance the contri- 
bution C\jr from the nonlinear term. Therefore, we need (1 — 2R 2 ) to 
be negative. This justifies the imposed lower bound R > l/y/2. 

4.2. Existence of a supersolution and a subsolution. Let e > 

be a small constant. Consider the function 

u+ = e{Ui + e 2 v a ) } 
where a is a constant to be chosen later. 

c 3r 3 Di(Ui + (Pv^D^Ux + £ 2 u )Ai(^i + e 2 v a ) 

{61) CU + — —E LV a — e — ■ = ryr . 

We have |-Djf a | < 2a and \DijV a \ < — , which combined with ( 12~T1) . give 
us 

V«+ < + -(1 - 2i? 2 ) + (d + £ 2 2a) 2 f ^ + e 2 ^ 

= I ( _ a (2i? 2 - 1) + Cf + 10ae 2 C 2 + 28a 2 £ 4 C 1 + 24a 3 e 6 

Here we used the notation Ci = K\ We replace a by a = /cCf, 
with a constant to be chosen later, to get 

(32) -£u + < C\ I - k(2R 2 - 1) + 1 



+ C\ (lOe^Cl + 28(e 2 kC 2 1 ) 2 + 24(e 2 kC 2 f) . 
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If e is small enough so that r\ := e 2 kC\ < 1, the equation ( 1321) becomes 



^£u+ < C\ (-^f- - k(2R 2 - 1) + 1 + 6277 



Hence, to have a supersolution, it suffices to find k and e that satisfy 
71 ■= £ 2 kC 2 < 1 and - k{2R 2 - 1) + 1 + 62r? < 0. 

If we take R Q so that R > R > l/v2, both inequalities are true for 
k = and e 2 C 2 < ^ = ^§^. Recall that C x = tfi||F|| C 4, 

therefore, for e and k so that 



j / 2/^2 j 2 

^llc«<^VHiT-- * = 2flFT md BSBo>1 

the function u + = e(U\ + e 2 kC\{r — — )) is a supersolution £u + < 
with boundary condition w + |an = eF. A similar argument shows that 
w_ = e{JJ\ — e 2 kC\{r — — )) is a subsolution with boundary condition 
U-\an = £F. If we denote by / = eF and ui = eU\, we just proved the 
following result. 

Lemma 11. Assume R > Ro > l/y/2 and denote by Q = R 2 \ B R , 
the plane with a hole of radius R at the origin. Let f be a function in 
C7 4 [0,27r] with the symmetries (J2J), U\ a solution to the linear equation 
Lu\ — 0, U\\dn — f and K% a constant (independent of f) so that 

\D Ul \ + \C\\D 2 Ul \ <#i||/||o*. 

Ife= II/IIc 4 [o,2tt] satisfies 



1 2R 2 - 1 



o<,= llfllc<^VT- 

then the function u + = U\ + Kfe 3 k(r — ^-), where k = 2 j^-i > satisfies 

£u + = g %3 (Du + )DijU + — £ ■ Du + + u + < in Q 
u + = f on dQ, 

and the function U- = U\ — Kfe 3 k(r — — ) satisfies 

Eu_ = g ij (Du-)D ij u- - £ ■ Du^ + u- > in Q 
u_ = f on dVt. 

From the bounds on the derivatives of U\ and the definition of u + 
and we have 
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Corollary 12. Let u + and w_ be defined as in the lemma above. Then 
there exists a constant K 2 depending only on ||/||c 4 an d Ro so that 

icr'Ki + \d u± \ + leipvi < k 2 , ten 

5. The Parabolic Equation d T u = £u 

In this section, we prove that there exists a solution to the parabolic 
equation d T u = £u for all time r G [0, oo) for well chosen initial and 
boundary conditions u and / respectively. 

To avoid confusion, let us fix notations: we denote by Di the or- 
dinary derivative with respect to z-th coordinate Xi or & of a spatial 
variable and by Dv the gradient of the function v with respect to spa- 
tial variables. The notation d T or d t is reserved for the derivative with 
respect to time. 

5.1. Initial condition and short time existence. For the initial 

condition to the parabolic equation, we choose a C 3 function u in fl 
that stays between the sub and supersolutions constructed in Lemma 
[Til and has bounded derivatives 

U- < Uq < u + , 

\\D$u \\ C 2(n) < Mo < oo- 

Moreover, we assume that uq satisfies the following symmetries and 
compatibility conditions 

(33) uo(r, 0) = -u (r, -9) = u {r, tt/N - 9), 

(34) u (R,9) = f(9), S(u )(R,9) = 0. 

Let us define the domains Ot = {x G R 2 | |x| > y/2(l-t) R} and 
Q T = {fat) | < t < T, \x\ > y/2(l -t) R}. We consider the 
following corresponding problem for the mean curvature flow 

(35a) d t v(x,t) = g ij (Dv)(x,t)D ijV (x,t), 

(35b) v(x, 0) = v (x),x G Oo 

(35c) v(x, t) = V2f( ; = =),t > 0, x G dO t . 

for some initial condition Vq(x). Standard parabolic theory assures the 
existence of a smooth short time solution v(x, t) to the boundary value 
problem (1331) with initial condition v (x) = v / 2w (-^). Moreover, if we 
denote the maximal time interval in which v exists by [0, to), we have 

(36) M t = sup \Dv\ + sup \D 2 v\ < oo, t G [0, t ). 
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For more details, we refer to [2j. 

In order to show that the solution v exists for all time t, we argue by 
contradiction and assume that t < oo then show that M t is uniformly 
bounded for all t G [0, t ). Standard parabolic theory implies that 
the solution v can then be continuated past t , which contradicts the 
maximality of [0, to)- The rest of this section is devoted to the proof of 
a uniform bound on M t . 

The function u defined by 

x 1 



(37) v(x,t) = y/2{l - t) „(__, -- ln(l - t)) 
satisfies 

(38a) d T u = g ij {Du)D ljU - £ • Du + u in Q x [0, r ) 

(38b) u(t 0) = Uo(0, e e n, r) = /onffix [0, r ), 

where we used the change of variables £ = , x , r = — Mn(l — £) 
and r = —\ ln(l - £ ). 

Lemma 13. Let u be a smooth solution of (138aj) in Q x [0, To) wi/t 
boundary condition (I38bl) . t/ien 

(39) max l-Dful + max I-DLmI < M(ri) if ti < r 

V ' nx[0,n] nx[0,n] K 



Proof. To each solution it to (j38|) corresponds a solution t> to the prob- 
lem ( 135]) via the formula fl37|) . This lemma is therefore a consequence 

of urn □ 



5.2. Bounds on the function u. Here, we show that a solution u to 
the problem (138]) stays between the subsolution w_ and the supersolu- 
tion u + : 

Lemma 14. Let u(£,t) be a C 2 solution to the problem f[3"8"j) . then 

(40) u _(0 < u(tr) < w + (0, > i?,0 < r < r Q . 

Proof. The proof is a slight modification of the proof of a Maximum 
Principle for parabolic equations. Let Ki be as in Corollary [12j For 
small e > 0, consider the function IV = ip(u- —u)—e, where if}(£, r) : = 
e -z\£\-( 1 + s ) r anc i 5 j s to be chosen later. Note that W becomes negative 
as |£| tends to infinity. We have d T ip = —(1 + 8)ip,Dii() = —e^ip, 



and Dijip 



■0. A computation shows that IV 
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satisfies 

(41) d T W < -8(W + e) - W ■ DW - e\£\(W + e) 
+ ^(g l >(Du_)-gV(Du))D ljU „ 

+ g ij {Du) (D^W + 2ej|r • DW + (W + e)^- 

Fix 7~i G [0, Tg). The estimates in Corollary [T21 and Lemma [T51 guarantee 
the existence of a large constant rj such that e~ e '^(|w_| + |w|) — e < — e/2 
for |£| >T],t& [0,ri). Denote by A R;V the annulus {i? < |£| < 77}. On 
the boundary of the cylinder A R ^ x [0, ri], we have 

W(£,r)<-e/2 |£| = r),r G [0,n] 

W(£,r)<-e |£| = R,r E [0, ri] 

W(f,0)<-e A^x{0}. 

Suppose that W > for some point in the cylinder Ar, v x [0, ri]. There 
is a first time r 2 G (0, Ti) for which there is a point £ with W(£, r 2 ) = 0. 
At this point (|,t 2 ), DW(|,t 2 ) = 0, -.D 2 W(£, r 2 ) is positive definite, 
and <9 r J¥(£,r 2 ) > 0. At (|,r 2 ), the equation dHJ yields 

< - £ 2 |^| + i){g ij {DuJ) - g ij (Du))D ijU ^ + eg ij (Du'" D "' 



The last term can be estimated by eg lj ^j^- < ^ + ^ + e 3 < As 2 . To 
bound the penultimate term, note that 



\d P J 3 ip)\ 



5ikPj + 5 jk Pi <ipiPjPk 



< 6 



1 + bl 2 (1 + M 2 ) 2 

therefore, |V>(^(£>U-) - 9 ij ( D u)) D ijU - \ < \2^[2e 2 K 2 at (£ 2 ,t 2 ) by 
the Mean Value Theorem. We obtain 

< -eS - e 2 \&\ + l2V2e 2 K 2 + Ae 2 < e{-8 + \2\/2eK 2 + Ae). 

The choice of S = 18eK 2 + Ae makes the last expression negative. 
This is not possible. We have therefore shown that W(!;, r) < for 
R < |£| < f]i T ^ [0) r i]- The argument stays valid for any rf > 77, thus 
r) - u(e, r) < ee^l+a+^^Jr r ) e n x [0 , r J. 

Letting e tend to 0, we obtain u_(£, r) < «(£, r), for (£, r) G Q x [0, ri]. 
Since T\ is arbitrary, we have 

U-(£,t) < u(£,r), for (£, r) G x [0,r ). 
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The inequality w(£, r) < «+(£, r) is proved in a similar way. □ 

5.3. Bound on the first derivative \Du\. We prove the bound in two 
steps. First, we use Lemma [T41 and interior estimates from the mean 
curvature flow to bound \Dv\. The interior estimate on \Dv\ gives us 
estimates for \Du\ for points staying farther away as time progresses. 
The second step is to consider an annulus Ar^ = {R < |£| < r/}, with rj 
large enough to have bounds for \Du\ on |£| = r\ from the first step. A 
maximum principle for the parabolic equation satisfied by \Du\ 2 yields 
estimates in the interior of the annulus, and therefore on the whole 
domain. 

Lemma 15. Let r) be a smooth solution to ( I38al) in the cylinder 
Q x [0, To) with boundary conditions f!38bj) . There is a constant Mi 



depending only on the boundary conditions Uq and f such that 
(42) \Du(tr)\<M 1 
for all times r for which u exists, and |£| > R 
Proof. Step 1. Away from the boundary. 

Let v(x, t) be the corresponding solution to the MCF defined by equa- 
tion fl37j) . We use the interior estimate below established by Ecker and 
Huisken: 

Theorem 16 (Theorem 2.3, p551 pQ). The gradient of the height func- 
tion v satisfies the estimate 



(43) y/l + \Dv(x ,t)\* < Ci(n) sup y/1 + \Dv(x, 0)| 2 

B p (x ) 

exp[C2(n)p~ 2 sup( sup v — v (xq, t)) 2 ] 

[0,T] B p {x )x[0,T] 

where < t < T, B p {xq) is a ball in Q and n is the dimension of the 
graph of v. 

The first factor is estimated by 

(44) \Dv \ = \Du \ < \\Du \\ c2m = M . 

To bound the second factor, we recall that «_(£) < u(£,t) < u+(£), 
therefore Corollary [121 gives us 

(45) \v(x,t)\ < K 2 \x\. 

The function v(x,t) is defined on the domain {(x,t) | < t < to, \x\ > 
a/2(1 — t) R}, in particular, it is defined for |x| > R\[2 at any time 
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< t < t . Let \x \ > 2R^/2 and p = \x Q \/4. The estimates (jSj) and 
(USD yield, 

y/l + \Dv{x ,t)\ 2 < C for \x Q \ > 2V2R and t G [0,t ). 
Note that Du(-yJL=, -\ ln(l - *)) = Dv(x,t). With the change of 
variables r = — ~ ln(l — t) and £ = , x = ^C, we get 

2 V ' S V^2(l-t) V2' to 

(46) |^(e,r)| < K 3 , for |f| > 2Re T , r G [0,r ). 
Withouth loss of generality, we can assume that K 3 > M = ||Z>Uo||c 2 (n)- 

Step 2: Up to boundary, using a Maximum Principle. 
A computation shows that the function w = \Du\ 2 satisfies the para- 
bolic equation 

d T w < g i3 (Du)D ijW + g^(Du)DijU D lW - f • Dw, 

where g ij {p) = <J y - jggjj and g% = ^g ij . Denote by Z { = g^(Du)D i:j u, 
then 

(47) d T w < g ij (Du)D ij w + Z l Diw - £ ■ Dw. 

Choose < T\ < tq. It follows from Lemma [13] that the coefficients in 
( 1471) are bounded on Q x [0, 7jJ. Moreover, from ( 1461) . we know that 
ty is bounded in the cylinder Q x [0, Ti]. A Maximum Principle for 
parabolic equations on unbounded domains (see Theorem 8.1.4 in [5] 
for example) implies that the function w in Q x [0, t\) is smaller than its 
supremum on the parabolic boundary [dVL x (0, Ti)] U [Cl x {0}]. Since 
T\ is arbitrary, the result is valid on the time interval [0, To) as well, 

sup w < sup w < max(K%, Mq). □ 

nx(o,r ) [dnx(o,T )]u[nx{o}] 

5.4. Bound on \D 2 u\ and \D 3 u\ away from boundary. 

Lemma 17. Let u be a solution of the problem (J3"8"j) inQx [0, t ), i/ten 
we /lave i/ie following bounds for the second and third derivatives of u: 

(48) \DIu(£,t)\ < e~ 5o C+ max |D| f u|, |£J > e(i? + 1), t G [0, t ) 

f2x[0,(5o] 

(49) \DIu(£,t)\ < e- 25o C+ max \DLu\, |f| > e(R+ l),r G [0,r ) 

nx[0,(5o] 

where 5q = min(T /2, 1) and C denotes different constants depending 
only on 5 , \\Du \\c2(n) and \\f\\c*. 
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Remark 18. The inequality ( HHj) does depend on r through the con- 
stant 5q; however the dependence is very loose. In the proof, we will see 
that Sq can be chosen to be min(l, s/2) for any time < s < r . There- 
fore, the estimate is uniform in time once we know that the solution 
exists past some small time s. 

Remark 19. maxQ X [ ,«5 ] \D^u\ and max^x^^j \D'^u\ depend on u on 
the right hand side, but we are primarily interested in finding a bound 
for r approaching r so there is no need for a better estimate near initial 
time. 

Proof. Let us tackle the second derivative first. The estimate for r < 5 
is immediate. For r > 5q, we use the change of variables t — 1 — e~ 2r 
and look at the function 

(J37D v(x, t) = V2(l-*) U (-^===, -\ ln(l - t)) , 

which satisfies the mean curvature flow (I35ap . We are interested in time 
t such that t > 1 - e~ 25 °. We can use the following interior estimate on 
the curvature from Ecker-Huisken [I] as long as v exists in the domains 
mentioned. Define t to be 1 — e~ 2l ~°. 

Lemma 20 (Corollary 3.2 [lj). Let < t x < t and rj > and < 9 < 
1. Fort G [0,£i], we have the estimate 



sup \A\ 2 (t)<c(n)(l-6 2 )- 2 [- + -) sup (l + \Dv\ 2 ) 2 

B ev (x ) \V tj B v (x )x[0,t] 



where B^xq) denotes the ball of radius rj centered at xq in the plane. 

First, we show that v exists in B v (xq) x [0, ti] for rj = 1, 6 = 1/2 and 
l^o I > V2(R + 1): if (x,t) G B v (x ) x [0,£i], then \x\ > \x \ — 77, and 
the corresponding f = . x satisfies If I > 'f '"" > v/2( / fl+1) ~ 1 > i2. 

Hence, for |a? | > y/2(R+ 1) and all t G [1 - e -M °,to), 

16 1 

|A|W)< sup |A| 2 (t)<L7-(l+ )L7<L7(^ ). 

-81/2(2:0) y 1 e 

Bounds on the second fundamental form A and on the first derivative 
Dv yield a bound on the second derivative of v: 

\D 2 xx v(x ,t)\<C(8 ), \x \>V2(R+l),te[l-e- 26o ,t ) } 

where C(5 ) denotes a different constant, also dependent on S . There- 
fore 

\D 2 r)| = >/2e"'-|DLi;(V2 e -^ J 1 - e^)\ 



(50) 



<e- r L7(<5 ), |e| >e T (R + l),re [5 ,r ). 
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In order to prove an estimate independent of time and outside a fixed 
annulus, note that for all T\ > 0, the function u(£, s) = u(£, s + Ti) is 
also a solution of d s u = £u. Moreover, \D^u(^,s)\ < Mi for (£, s) G 
Q x [0, r — Ti) so we have an estimate similar to (150]) in this case also 

s)| < e~ s C(5 ), |e| > e s (i? + 1), s G [<5 , r - n). 

Let s = 5q and r = n + s, the inequality above yields 

|^>(e,r)| < e- 5o C(5 ), |f| > e 5o (i? + 1) 

for all 7~i > 0, i.e. for all r > 5 . 

The bound on the third derivative is proved in a similar fashion using 

/ 1 l x 2 
sup \VA\ 2 (t) < c — + - 

Bfl„(a:o) V 7 ? * 

with the constant c depending on 9 and sup B ( xo ) x [o,t](l + l-^'l 2 ) from 
Theorem 3.4 in pQ. □ 

5.5. Holder continuity of D^u near boundary. Let the time T\ < 
t and denote by 

A R>Rl the annulus {£R 2 | R < |£| < Ri}, 

QR,Ri,n the domain A R)Rl x (0,Ti), 

T RiRljT1 the set (A flijRl x {0}) U ((dB Rl U x [0,n]), 

Mi = max iZVul 
nx[o,T ) 

For a domain Q G R 2 x [0, oo), denote by C 2,1 (Q) the set of all contin- 
uous functions in Q having continuous derivatives D^, D 2 ^u,d T u in Q 
and by C a,a l 2 (Q) the set of functions on Q that are a-Holder continu- 
ous in £ and a/2- Holder continuous in r. 

We use the following theorem by Ladyzenskaja et al. 

Theorem 21 (Theorem 2.3, p533 [6]). Let u be a solution of equation 

(51) d T u — ay(£, r, w, D^u)D^u + a(£, r, w, Z^w) = 

belonging to C 2,1 (Q R)RltT1 ) . Suppose that equation (J5TI) parabolic with 
respect to u(£,t), i.e the inequality 

V\C\ 2 < tkj(.Z,TMt,T),D{u(Z,T))CiCj < HCI 2 , v,»>o 

is fulfilled and the functions r, u, D^u) and a(£, r, u, D^u) are con- 
tinuous and continuously differentiate with respect to all of their ar- 
guments in the domain 

{{Z,t) G Q r ,r 1iT i, \u\ < max |u(£,t)|, \p\ < M x = max |w € (£,r)|l 
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Suppose also that there is a constant M 2 so that 



max max 



aij(£,T,u(£,r), Dgi(£,r))\; \a\ 



dai 



da 



dr 



da,; 



dr 



da; 



da 



dD^u 
da. 



da 



du 



da 



du 



dD^u 
<M 2 , 



is estimated 



Then for some a > 0, the quantity \\D£u\\ca,>x/2(Q RR ) 
by a constant depending only on Mi, the quantities rj,fj, and M2, the 
norm IM|c 2 . 1 (r R .R 1 n ) an d the norm in C 2 of the functions defining the 
boundary of Ar^. These same quantities also determine the exponent 
a. 

The coefficients of the parabolic equation 

( Egg) d T u = g ij (Du)D ij u-^- Du + u 

m Qr,Ri,ti satisfy the hypotheses of Theorem |2"T1 We therefore have a 
bound on \\D^u\\ ca ,a/2^ RiRliTl x[o,Ti]) tnat depends on \\f\\c*, R, Mi and 
||'u||c a ' 1 (ri i i{ 1 T )• From Lemma [T71 we can estimate this last quantity 
uniformly in time, therefore, 

(52) \\ d zu\\c<*><*/2(a r , Ri x[o,to)) <K 4 <oo, 
where only depends on time through 5q = min(l, tq/ '2). 

5.6. Holder continuity of D^u near boundary. The coefficients 
g^(Du) of (I38al) are in C a > a / 2 (Ar^ x [0,r )) hence we can apply 
standard theory for linear parabolic equations with Holder coefficients 
(Theorem 5.2 p320 [6] for example) to obtain 

(53) IHIc^.a/a+i^^x^ro)) < -^5 < °0- 

The Remark [T5] explains the loose dependence of our bounds on S and 
t . Let us emphasize it again: if the solution u is known to exist past 
a small time s > 0, we can choose Sq = s/2 and the constants in all the 
estimates following (HHl) do not depend on r > s. 



5.7. Continuation of the solution. Combining 
have 



ID and flM]), we 



|-D^w||c«,«/2(nx[o,r )) 



~T" J J d T 1l 1 1 (ja,a I 



< K R < 00. 



*/ 2 (Ox[0,t )) 

Moreover, Lemma [TBI gives us ||-D§w|| < Mi < 00 on Cl x [0,to). These 
estimates can be transformed into uniform bounds on the first two 
derivatives of the corresponding solution v to the mean curvature flow 
( I35al) for time t G [0,t )- Therefore the solution v has a continuation 
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existing for time t G [0, £3) with £3 > to- This contradicts the maxi- 
mality of to, and proves that the solution u to (1551) exists for all time 
r G [0, 00). Since none of the bounds on u, Du or D 2 u depended on 
t , we have 

(54) «_(£) < «(£, r) < (£, r) G x [0, 00), 

(55) |£) c u(e,r)|<Mi, (£, r) G O x [0, 00), 

(56) ||-D^w|| cq W2(q x [ 0)Oo)) + ||9rM||c Q ^/ 2 (Qx[0,oo)) < ^6- 



6. A Solution to the Elliptic Equation Su=0 

Here, we show that there is a subsequence of u(-, r n ) that tends to a 
solution to the elliptic equation £u = as r n — ► 00. 

Let us denote by M T the graph of u(-, r) over fi, by X = (£, r)) 
the position vector of a point on M T . From the equation d T u = £{u), 
the variation of X with respect to time is 

d T X = (H + X ■ u)y/l + \Du\ 2 e 3 =: Z, 

where e 3 is the unit vector (0, 0, 1) in R 3 . It will be useful for what fol- 
lows to decompose e3 and Z into their tangent and normal components 
with respect to M T : 

v 

e 3 = v + 



1 + \Du\ 2 ' 



Z = Z T + Z 1 - = (H + X ■ v)y/l + \Du\ 2 v + (H + X ■ v)v 

where v = ^^^ {D 1 u,D 2 u, \Du\ 2 ). 

Consider the functional J(u) = j M e~^ 2 ^ 2 dH, where dH is the 
Hausdorff measure on M T . We have 

= f UdwMZ 1 + div M Z T ) - X ■ d T x] e-l*! 2 / 2 dH 

= / e~W 2 / 2 (-H(H + X-v)) dH+ I e~\ x \ 2 l 2 div M Z T dH 

J M T J Mr 

-fx- d T Xe~\ x \ 2 / 2 dH 

J Mr 

= - [ e-\ x ?l 2 {H + X-v) 2 dH+ [ e-\ x \ 2 ' 2 Z T ■ rj ds, 

J Mr JdMr 
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where 77 is the outward co normal unit vector to dM T . On the circle 
r = R, = d T u = (H + X ■ v) y/l + \Du\ 2 therefore 

e -\x?/^ H + X -v) 2 dH<0. 

Writing in terms of d T u and as an integral over f2, we obtain 

dr J n ^\ + \Du\ 2 

By definition, J(u(-, r)) > for all time r G [0, 00), so 

f°° f (d T u) 2 i«i 2 + . 2 
Jo 7n ^/l + \Du\ 2 

and 

(57) hm / / — ===e ~ d£ = 0. 

n ^°°Jn-i Jn a/1 + \Du\ 2 

Define w n (£, r) = w(£,?2 + r) for (£, r) G x [—1,1]. For any con- 
stant p > R, the uniformly bounded in C a+2 ' a/2+1 ((B p HQ) x 
[—1,1]). We can therefore extract a subsequence that converges in 
C0+ 2 '0/ 2+1 ((5 p nfi) x [-1,1]), < /3 < I. A Cantor diagonal argu- 
ment gives us a further subsequence u n% and a function « M such that 
— > Moo as i — > 00 in the norm C /3+2 ' /3 / 2+1 ((£? nfi) x [—1,1]) for ev- 
ery 00 > p > R. Since every u n satisfies d r u n = £u n in ft x (— 1, 1) and 
w n(£> t) = / on <9-Br x [—1, 1], the limit function Uoo inherits the same 
properties. On the one hand, equation (1571) and the uniform bounds 
( |54|) and (1551) tell us that d T u n ^ — *> pointwise. On the other hand, we 
know that d T u ni — *> d T Uoo pointwise. Hence, d T Uoo = and -u^ satisfies 
Euoq = in Q and = /, dQ. We have therefore proved 

Lemma 22. Let 2 > R > ^ and N > 5. Define Q = {£ G R 2 | |£| > 

i?} to 6e i/ie p/ane a /10/e of radius R. There is a E\ > depending 
on R and N such that for any f : [0, 2n] — > R || /" || c 4 [o,2tt] = £ < £1 
and satisfying the symmetries f{6) = —f(—0) = f{it/N — 6), there 
exists a smooth solution u to the Dirichlet problem 

© £{u) = g ij (Du(0)D ijU (£) - £ ■ Du(£) + u(£) = 0, £ G fi, 
(SD u = f on dB R . 

In addition, We can choose the constant e± = Ei(Rq) uniformly for all 
R such that 2 > R > R > 1 /y/2. 
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7. Uniqueness 

For the sake of clarity let us assume R > a/3/2 in this section. The 
same reasoning; also works in the more general case R > 1/ a/2 although 
the estimates are more involved. 

Let / be our boundary condition, and suppose that ||/||c 4 (9c) < £ i, 
where E\ is given by Lemma [221 

Theorem 23. There is an e 2 < £\ so that, if ||/||c 4 < £ 2 and if u and 
u are two solutions to 

(58a) Su = g ij (Du)D lj u -£-Du + u = 0inQ 

(58b) u = f on dn 

with U- < u,u < u + then u = u. 

Before we start the proof, we need some a priori estimates on the 
solutions u and u. 



7.1. A priori estimates. Because the solutions u and u above are not 
dependent on time, they satisfy the parabolic equation d T u = £u as 
well. Estimates similar to the ones from Lemmas [TS] and [T7] therefore 
hold in this case also, with some small modifications. If u is a solution 
to the elliptic equation f)58ap . then the function 

is a solution to the mean curvature flow on Qi = {(x,i) | < t < 
l,\x\> y/2(l -t)R}. The estimate 

(USD y/l + \Dv(x ,t) | 2 < Ci(n) sup a/1 + \Dv(x, 0)| 2 

B P {x ) 

exp [C 2 (n)p~ 2 sup ( sup v — v (x , t)) 2 ] 

[0,T] B p (x )x[0,T] 

where < t < T,B p (xo) is a ball in Q and n is the dimension of 
the graph of v , is valid. The last factor of the right hand side can 
be bounded as in the proof of Lemma [151 however, we do not have 
an initial condition independent of v in this case, so right now, our 
estimates depend on the value of Dv(-, 0) over a bounded set. Taking 
\xq\ = |a/2_R and p = |xq|/3, we obtain 

\Dv(x ,t)\<C sup a/1 + \Dv{x, 0)| 2 , 1 > t > 0, 

2 v / 2i?>|x|> v / 2i? 
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hence, by the change of variable x = £y2(l — t), 

\Du{£)\<C sup y/l + \Du(()\ 2 , \£\>R 

2R>\C\>R 

An argument similar to the one in the proof of Lemma [17] gives us 
bounds on higher derivatives of u away from the boundary 

\D 2 a u(0\<C, \£\>e(R + l), 

\D^u(0\<C, \£\>e(R+l), 

with the constant C depending on sup 2 R>^|>/j y/l + \Du(()\ 2 . Stan- 
dard elliptic theory on quasilinear equations (see Theorem 13.4 or 13.7 
[3] for example) assures the existence of constants C and a > de- 
pending only on /, R and sup 2i?> |£| >J? a/1 + \Du(()\ 2 such that 

\\ Du \\c a (A R , 2R ) < C, 

where A R2 r is the annulus {R < |£| < 2R}. The C a Holder norm of 
the coefficient g % i(Du), seen as a function of £, is uniformly bounded 
over the whole domain fl. We have a bound on ||^|^w||c»(Ab 2 r) ^y 
elliptic theory; therefore ||f |^w||c a (n) is also bounded. 

7.1.1. The gradient Du achieves its maximum on the boundary dfl. 
Equation (1471) from section implies that w = \Du\ 2 satisfies 

(59) < g ij (Du)D ij w + Z k D k w - £ ■ Dw =: Cw 

where Z k (x) = g k J (Du)DijU and g k {p) = ^f-(p)- We will derive the 
following theorem as an application of the maximum principle: 

Theorem 24. Assume u G C 2 {Vl) fl C°{pl) is a solution of (I58I) . then 
\Du\ achieves its maximum on the circle r = R. 

From the discussion above, \Z\ is uniformly bounded and L is a 
uniformly elliptic operator on Q. Moreover, the coefficients of L are 
Holder continuous with bounded Holder norm on any compact set of 

a 

Lemma 25. If\Wh\ < b for some constant b, there is a function ip(r) > 
for which C{ip) < and ip(r) — > oo as r — ► oo. 

Proof. We consider positive increasing functions ip(r). The operator 
C(ip) is estimated by 

= g ij (Du) ((^ - + + U'Z -i-r^' 

<^(Du)^"+(l-r+\Z\^'. 
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We have uniform bounds on \Du\ and \Z\, therefore 

CQ>) < (g^(Du)^f) (f + ^(i - r + C){1 + b 2 )). 

The solution to the differential equation ip"+ip' (- + C — r) (l+b 2 ) = 
is given by 

ln(^'(r)) - Hi>'(R)) = (1 + b 2 ) (- ln(^) - C(r - R) + -{r 2 - R 2 ] 

\ R 2 



1 J>(r)-1>(R)=iJ/(R) I i-) 1+b e (-c(s-R)+Us*-R*Mi+b*) d 



s 



s. 



For the choice of boundary conditions if>(R) = 4>'(R) = e > 0, the 
function ip and its derivative ifi' are positive for r > R. The last thing 
to verify is that if) tends to infinity as r approaches infinity. This is 
true since for large r, the dominant term in the expression for if}' is of 
order e r . □ 

Proof of TheoremWA Let us consider the function w — rjip, with ip as 
in Lemma [231 and w = \Dh\ 2 . This function is negative for large values 
of r. From equation (|59l) . we have 

C(w — rjip) > 0, 

for any positive constant rj. By the maximum principle, w—rjip does not 
achieve an interior maximum. Since the function w — rjip goes to — oo 
as r goes to infinity, it has to achieve its maximum on the boundary 

w {0 — W&dCl) — maxw — t]e, £ efl,i)>0 

Letting r\ tend to 0, we obtain u>(£) < max^nw for £ E Q. □ 

Our subsolution and supersolution provide barriers for the function 
u, therefore 

|£H<C||/|| C4 

with the constant C independent of / and u. With this bound, we can 
run through the arguments in the beginning of this section again and 
obtain estimates for the second derivative D 2 u that are independent of 
the solution u. 

7.2. Proof of uniqueness. 

Proof of Theorem\^ Let u and u be two solutions as described in the 
hypotheses. We substract equation (158aj) for u from the equation for u 
and use the notation w = u — u to obtain 

g lj (Du)D ijW + [g ij (Du) - g ij (Du)]D i:j u - £ ■ Dw + w = 0. 
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The Mean Value Theorem implies that 
g ij {Du)D ij w + [gi j (Du + 9(Du - Du))}D i:j uD k w - £ ■ Dw + w = 0, 

for some 6 6 (0, 1) and where g^(p) = for p e R 2 . Denoting by a 

the vector with components = [g^ (Du + 6(Du — Du))]DijU, k = 1, 2, 
we have 

(60) C(w) = g ij (Du)D ljW + a ■ Dw - £ ■ Dw + w = 0. 

A computation shows that g^(p) < Q\p\ if \p\ < 1. From Section ITTLj we 
know that there exists a constant C so that \Dijii\ < C in Q. Moreover, 
the maximum of the first derivative of any solution to (158aj) is achieved 
on the boundary dQ hence max^ \Du\ and maxQ \Du\ are controlled by 
e = II /lie 4 an d go to zero as e tends to 0. The same fact is true for 
rj = sup \a\ by the definition of a. 

The maximum principle does not apply immediately to flBUl) because 
of the positive coefficient in front of w. To circumvent this, we use 
a positive supersolution ip to (!60j) which grows faster than linearly at 
infinity. The existence of such a supersolution ip is given below and is 
the main difficulty in this proof. 

Assuming such a positive supersolution ip exists, we define w = w/ip 
and look at 

C(w) := C(ipw) = Cw = 0. 

The coefficient in front of w in the explicit expression of C(w) is now 
C(ip), which is non positive, so the maximum principle applies to C. It 
implies that for any sub domain Q' of Q, 

sup w < sup w. 

From the hypotheses, w = on dQ and from the growth of ip at infinity, 
|u>(£)| tends to as |£| tends to infinity. Choosing the domains Q' to be 
increasingly large annuli, we obtain that sup^w) = 0. Since w = w/ip 
and ip > 0, sup n u> = 0, therefore u < u in the entire domain Q. 
Switching the roles of u and u we show that u < u and conclude that 
u = u in Q. 

We are left to prove the existence of the function ip. 

Claim 26. For R > V3/2 and 82 small enough, there exists a function 
ip on Q satisfying 

l£l 

ib > 0, Cib < 0, and lim = 0. 
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The supersolution is found among functions depending on the radius 
r only, ip = i/>(r), that are increasing ip'(r) > 0. The operator £ applied 
to a radial function is estimated by 



C^)=g»{Du) ( 



r + 77 

r 



Let us consider functions of the form 

iP( r ) = r a - br a ~ 2 

with constants 2 > a > 1 to be chosen later and 6 = 3/4. This function 
is positive on the set r > R provided b < R 2 . To simplify the notations, 
we work with the case 6 = 3/4 and imposed R 2 > 3/4; however, the 
following argument is also valid in general for R 2 > b > 1/2 with 
slightly more subtle estimates. The first and second derivatives are 
given by ip' = ar a ~ l — b(a — 2)r a ~ 3 > and ip" = a (a — l)r a ~ 2 — b(a — 
2) (a - 3)r a " 4 . We then have 



g ij (Du)^-ip" <a(a-l)r a - 2 



1 + \Du\ 2 



b(a - 2) (a - 3)r 



a-4 



By choosing e small enough, we can assume that \Du\ 2 < 1/8. A 
computation gives 



„a-2 



£(V0 < (1 - a)r a + riar a ~ Y + (a 2 + ba- 36)r 

+ 776(2 - a)r a ~ 3 + 6(2 - a)(l - -(3 - a))r a ~ 4 

9 

The coefficient of r a ~ 2 is negative ifl<a<9/8< ~ 3+v ^ anc j 
6 = 3/4. Such a choice of a also guarantees that the coefficient of r a ~ A 
is negative. Therefore 

£(V0 < r a " 4 (^(l -a)r 4 + 77ar 3 + 77^(2 -a)r) 
(61) < r a ~ 3 ((l -ct)r + 2r/) 

since |r < r 3 . The right hand side of (ED) is negative if we choose e 
small enough and a so that 
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8. Remarks 

Let || /||c 4 < e an d let u be the unique solution to (1581) such that 
U- <u<u + . 

Remark 27. Because of uniqueness, the function u has to satisfy the 
symmetries 

p u(r, 6) = -u{r, -6) = u{r, n/N — 6) for r > R,9 E [0, 2tt). 

Taking e — £ 2, Theorem [1] from the introduction is then a corollary of 
Lemma [221 and Theorem [231 

Remark 28. From the explicit formulas in Lemma [TTl we have 

U- ~ Mi — £ 3 fc < U < M + ~ Mi + £ 3 fc, 

where the constant C and the function vc do not depend on e. Recall 
that U\ is the constructed solution to the linear problem Cu\ = 0, 
= f. This implies that D T u = D r u\ + 0(e 3 ) on the boundary 
circle r = R, in other words, the radial derivative of u on the circle 
r = R is the radial derivative of the solution to the linear problem with 
an error of order e 3 . 

Remark 29. Considering equation ([50]) with |£| = e T (R + 1) gives us 
l^>(e)l<||, \Z\>e(R+l). 

Since \Du\ is uniformly bounded in Q, the above estimate implies a 
bound on the mean curvature of the graph, and therefore on X ■ v: 

\X-v\ = \H\<^ \Z\>e(R+l). 

The graph of the solution u is therefore asymptotic to a cone at infinity. 
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